Late Decaying Dark Matter, Bulk Viscosity and the Cosmic Acceleration 
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We discuss a cosmology in which cold dark matter begins to decay into relativistic particles at 
a recent epoch (z < 1). We show that the large entropy production and associated bulk viscosity 
from such decays leads to an accelerating cosmology as required by observations. We investigate the 
effects of decaying cold dark matter in a A = 0, flat, initially matter dominated cosmology. We show 
that this model satisfies the cosmological constraint from the redshift-distance relation for type la 
supernovae. The age in such models is also consistent with the constraints from the oldest stars and 
globular clusters. Possible candidates for this late decaying dark matter are suggested along with 
additional observational tests of this cosmological paradigm. 

PACS numbers: 98.80.Es, 95.35.+d, 95.30.Cq 



I. INTRODUCTION 

Understanding the nature and origin of both the dark 
energy [l| and cold dark matter Q constitutes a signifi- 
cant challenge to modern cosmology. The simplest parti- 
cle physics explanation for the cold dark matter is, per- 
haps, that of the lightest supersymmetric particle, an ax- 
ion, or a heavy (e.g. "sterile") neutrino. The dark energy, 
on the other hand is generally attributed to a cosmolog- 
ical constant, or possibly a vacuum energy in the form 
of a "quintessence" scalar field 0, [1] or fc-essence 0, Q 
which must be slowly evolving along an effective poten- 
tial. See for a review. In addition to these explana- 
tions, however, the simple coincidence that both of these 
unknown entities currently contribute comparable mass 
energy toward the closure of the universe begs the ques- 
tion as to whether they could be different manifestations 
of the same physical phenomenon. Indeed, suggestions 
along this line have been made by many [3, [^, [l^, UH • 
See |12| for a recent review. 

In Ref. [Tsj yet another mechanism was considered by 
which a dark-matter particle could produce the cosmic 
acceleration. In that work it was shown that entropy 
production and an associated bulk viscosity could result 
from a decaying dark-matter particle. Moreover, that 
bulk viscosity would act as a negative pressure similar 
to a cosmological constant or quintessence. However, in 
that paper it was shown the decay alone was not suf- 
ficient to produce the observed cosmic acceleration. It 
was proposed, however, [l^ that some, but not all of 
the desired cosmic acceleration could be accounted for if 
the particle decay was delayed by proceeding thorough a 
cascade of long lived intermediate states before the final 
entropy-producing decay. 

In this paper we expand on the hypothesis that the 
dark energy could be produced from a delayed decaying 
dark-matter particle. Here, we show that a dark-matter 
particle which, though initially stable, begins to decay to 
relativistic particles near the time of the present epoch 
will produce a cosmology consistent with the observed 



cosmic acceleration deduced from the type la supernova 
distance-redshift relation without the need for a cosmo- 
logical constant. Hence, this paradigm has the possibility 
to account for the apparent dark energy without the well 
known fine tuning and smallness problems |12i] associated 
with a cosmological constant. Also, for the model pro- 
posed herein, the apparent acceleration is a temporary 
phenomenon. This avoids the difficulties fT3| in accom- 
modating a cosmological constant in string theory. 

The idea of delayed dark matter decay is not new. It 
was previously introduced [isj as a means to provide an 
Om = 0.1 — 0.3 without curvature ( fltot = 1) by hiding 
matter in weakly interacting relativistic particles. Here, 
we point out that such a cosmology not only allows for a 
flat cosmology with low apparent cold dark matter mat- 
ter content, but can produce an accelerating cosmology 
consistent with observations. We show that the bulk vis- 
cosity produced during the decay will briefly accelerate 
the cosmic expansion as matter is being converted from 
nonrelativistic to relativistic particles. This model thus 
shifts the dilemma in modern cosmology from that of 
explaining dark energy to one of explaining how an oth- 
erwise stable heavy particle might begin to decay at a 
late epoch. 

In the next section we summarize the cosmology of late 
decaying dark matter and its associated bulk viscosity. 
In the following section we discuss candidate particles 
for such decays, and in Section IV we present fits to the 
supernova magnitude-redshift relation which show that 
these data can be reasonably well fit in a flat A: = 0, A = 
cosmology with recent dark-matter decay. We summarize 
the constraints that the supernova data places on the 
properties of the decay along with independent constraint 
from the ages of oldest stars and globular clusters. 



II. COSMOLOGY WITH BULK VISCOSITY 

Numerous papers have appeared in recent years which 
deal with the subject of Bulk viscosity as a dark energy 
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[1, [l^. What is needed, however, is a physical model for 
the generation of the bulk viscosity. Below, we consider 
one possible means to produce a bulk viscosity in the 
cosmic fluid by decaying dark matter particles. To begin 
with, however, we first examine the effects on the cosmic 
acceleration of the bulk viscosity. For this purpose we 
utilize a flat (fc = 0, A = 0) cosmology in a comoving 
Friedman- Robertson- Walker metric. 



dr^ + r^dO^ + sin^ 



(1) 

for which U° = l,W = 0, and U)^ = 3d/ a. 

We consider a fluid with total mass-energy density p 
given by. 



p = PDM + Pb + Ph + Pi + pi 



(2) 



where pb is the baryon density, while pdm is the contri- 
bution from stable dark matter, p^ is the energy density 
in the usual stable relativistic particles (i.e. photons, 
neutrinos, etc.), and we denote the relativistic particles 
specifically produced by decaying dark matter as pi (al- 
though these products may well be normal neutrinos). 

Prior to decay, by fiat pi = pi{Q) = and the other 
terms in Eq. ^ obey the usual relations as given by the 
the conservation condition T^!'^ = and their respective 
equations of state, i.e.. 



dpi d 
— = -3- +k) 
dt a 



(3) 



where Pi the partial pressure from each species, so that 



Pdm 



PDAliO) 



Ph 



Ph{0) 



Pi 



PiiO) 



(4) 



Thus, for large C it is possible for the negative pressure 
term to dominate and an accelerating cosmology to en- 
sue. It is necessary, therefore, to clearly define the bulk 
viscosity for the system of interest. 

The Friedmann equation derives from the pv — 00 
component of the Einstein equation. Therefore, it does 
not depend upon the effective pressure and is exactly the 
same as for a non dissipative cosmology, i.e. 



H'' = - = ^-^Gp 



(10) 



where p represents to total mass-energy density from 
matter and relativistic particles (Eq. [2]). Even so, the 
bulk viscosity from particle decays can briefly affect the 
cosmic acceleration by producing a temporary condition 
of nearly constant p. 

Once decay begins, p remains as given above in Eq. ([2]). 
However, the conservation equations give new equations 
for energy densities in decaying matter and produced rel- 
ativistic particles. For pi and ph we then have. 



dph 
dt 



-^-{Ph +Peff) ~ ^Ph 



and 



d'Pi a"- I . \ . \ 



(11) 



(12) 



Denoting td as the time at which decays begin, for t > td 
we have the following analytic solutions: 



Ph = ^Ph(irf)e-(*-*-'/^ 



(13) 



where pi{0) denotes the initial energy densities at some 
arbitrary start time. 

We begin our models well into the time of the matter 
dominated epoch. Hence, p^iO) — aT^ w and the 
cosmology is nearly pressureless. However, once decay 
begins, the total energy density in relativistic particles 
P-y + Pi is not negligible even at the present epoch, neither 
is the pressure. Thus, we have: 



P = ^[P7 



Pi\ 



(5) 



With the introduction of bulk viscosity produced (as 
described below) from decaying dark- matter particles, 
the energy momentum tensor then becomes [13] 



Too 
T ■ 



p 





P 



3C^)5., 



(6) 
(7) 

(8) 



where this last equation shows that the effect of bulk 
viscosity is to replace the fluid pressure with an effective 
pressure given by. 



Pcff = P - C3- 



(9) 



and 



Pi{td) 



Ph{td) 



<t'-i^)/-a{t')dt' 



Pbv 



(14) 

where pBV is an effective dissipated energy in light rela- 
tivistic species due to the cosmic bulk viscosity. 



Pbv 



9 / at') 



a{t'Ydt' 



(15) 



The total density for the Friedmann equation will then 
include not only terms from heavy and light dark mat- 
ter, but a dissipated energy density in bulk viscosity. The 
introduction of this pBV term can lead to a cosmic ac- 
celeration as we shall see, but first, for completeness, we 
summarize the derivation of the bulk viscosity coefficient 
C given in [l^ and show how it results from the delayed 
decay of interest here. 



A. Bulk Viscosity Coefficient 

Bulk viscosity can be thought of [13, [H, [H, ^ as 

a relaxation phenomenon. It derives from the fact that 
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the fluid requires time to restore its equilibrium pres- 
sure from a departure which occurs during expansion. 
The viscosity coefficient ( depends upon the difference 
between the pressure p of a fluid being compressed or 
expanded and the pressure p of a constant volume sys- 
tem in equilibrium. Of the several formulations Il9| the 
basic non-equilibrium method ^231 is identical |l8l| with 
Eq. dH). 



C3--Ap , 
a 



(16) 



where Ap = p — p is the difference between the constant 
volume equilibrium pressure and the actual fluid pres- 
sure. 

In Ref. [18] the bulk viscosity coefficient was derived 
for a gas in thermodynamic equilibrium at a temperature 
Tm into which radiation is injected with a temperature 
T and a mean pressure equilibration time Tg. The solu- 
tion for the relativistic transport equation ^2l'| can then 
be used to infer [l^ the bulk viscosity coefficient. More 
specifically, the form of the pressure deficit and associ- 
ated bulk viscosity can be deduced from Eq. (2.31) of 
Ref. [l^ which we have generalized slightly and write as, 



Ap ■ 



T) 



Ap. 



7 ' c 



(Zdp\ 
\dp) 



where the subscript n denotes a partial derivative at 
fixed comoving number density. The factor of 4 on the 
r.h.s. comes from the derivative of the radiation pres- 
sure p ~ T'' of the injected relativistic particles, and 
the term in square brackets derives from a detailed solu- 
tion to the linearized relativistic transport equation [2ll | . 
This term guarantees that no bulk viscosity can exist for 
a completely relativistic gas (for which dp/dp = 1/3). 
In the cosmic fluid, however, we must consider a total 
mass-energy density p given by both nonrelativistic and 
relativistic components. 



B. Pressure Equilibration Time 

The timescale Te to obtain pressure equilibrium in an 
expanding cosmology from an initial pressure deficit of 
Ap(0) can be determined from. 



Ap(t) 
Ap(0)" 



dt 



(18) 



As in Ref. [l8[, in the present context we also have a 
thermalized gas into which relativistic particles at some 
effective temperature are injected. There are, however, 
some differences. For one, Eq. ([T7)) was derived under the 
assumption of a short relaxation time Te so that only the 
terms of linear order in were retained in the solution 
to the transport equation. In what follows we will keep 
this form of the solution even for long relaxation times 
with the caveat that this is only deduced from a leading 



order approximation to the full transport solution. We 
will, however, also consider a phenomenological analysis 
of the effects of terms of higher order in Te . 

Another difference in the present approach involves the 
nature of the pressure equilibration time Te. Indeed, this 
term contains the essential physics of the bulk viscosity. 
There are in principle two ways in which pressure equi- 
librium can be restored. One is from particle collisions 
and the other is for simply all of the unstable particles to 
decay. That is, we can write an instantaneous pressure 
restoration lifetime t as, 



1 



1 



'^decay '^coll 



(19) 



For the cosmological application of interest here the 
mean collision time TcoU — l/{nac) for weakly interacting 
(or electromagnetic) particles is very long (many Hubble 
times) and can be ignored. Hence, one only need consider 
the timescale to restore pressure equilibrium from the 
decay of unstable nonrelativistic dark matter. That is, at 
any time in the cosmic expansion the pressure deficit will 
be 1/3 of the remaining mass-energy density of unstable 
heavy particles. Hence, we replace p^/S with the pressure 
deficit from remaining (undecayed) mass energy ph/3 in 
Eq. (fTT)) and write. 



Ap ■ 



1 



,dp 



ax" 



(20) 



A form for the equilibration time Tc in the expanding 
cosmology can then be obtained from Eq. (fTS]) by setting 
T = Tdecay, and approximating H ^ d/a ^ constant. 
This gives, 



(21) 



[l + 3(a/a)T] 



Note, that the factor in the denominator acts as a limiter 
to prevent unrealistically large bulk viscosity in the limit 
of a large t. 

Following the derivation in 18], and inserting Eq. ([^U]) 
in place of Eq. pT|) , the form for the bulk viscosity of the 
cosmic fluid due to particle decay was deduced in [l^ to 
be, 



c 



1 - 



Pi +P7 



(22) 



where the square of the term in brackets comes from in- 
serting Eq. ([20]) into the linearized relativistic transport 
equation of Ref. [21I] . Note that equation implies a 
non- vanishing bulk viscosity even in the limit of long (in- 
flnite) relaxation times t as long as the total mass energy 
density p in the denominator is comprised of a mixture of 
relativistic and nonrelativistic particles so that the term 
in square brackets does not vanish. This apparent con- 
tradiction arises from keeping only the linear terms in the 
transport equation. Hence, one should be cautious about 
using this linearized approximation in the long pressure 
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relaxation-time limit. Even so, a more general derivation 
has been made [22| which shows that, even in the limit of 
interest here of a long radiation equilibration time there 
is a non-vanishing bulk viscosity consistent with experi- 
mental determinations. As an indicator of possible higher 
order effects in the relaxation time, we consider replacing 
Te in Eq. [22] with 



(23) 



where a or C is a parameter to be adjusted to fit the 
cosmological data. 



C. Evolution of Energy Densities 

To illustrate how an accelerating cosmology arises we 
plot the evolution of various energy densities in Figure 
[T] for a model with r = 20 Gyr and the onset of dark- 
matter decay at Zd = 0.3 {td ~ 8 Gyr). In this figure 
the dashed line shows the evolution of the energy den- 
sity in relativistic particles denoted pr = p-y + pi- The 
dot-dashed line shows the evolution of the total matter 
density pM ~ Pdm + Pb + Ph- From this figure it is 
evident that once the late decay begins, the sudden in- 
crease in radiation plus dissipative bulk viscosity leads to 
a finite period of nearly constant total mass-energy den- 
sity from the onset until well past the present epoch (at 
t = 11 Gyr). This mimics a A-dominated {p ~ constant) 
cosmology until nearly all of the unstable nonrelativis- 
tic dark matter has been converted to radiation. Af- 
terward, the total energy density continues to diminish 
as a radiation-dominated gas and eventually becomes a 
simple flat cosmology dominated by the remaining stable 
dark matter and baryons. 



III. CANDIDATES FOR LATE DECAYING 
DARK MATTER 

To avoid observational constraints, the decay products 
of any cold dark-matter particle must not be in the form 
of observable photons or charged particles. Otherwise 
the implied background in energetic photons would have 
been easily detectable [l^ . Neutrinos or some other light 
weakly interacting particle would thus be the most natu- 
ral products from such decay. We now summarize some 
suggestions for this type of decaying particle. 

A good candidate considered in [131 is that of a heavy 
right-handed (sterile) neutrino. Such particles could de- 
cay into into light i^g, v^, v-r "active" neutrinos ,24;] • In 
this case, the limitation of decays to non-detectable neu- 
trinos places some constraints on the sterile neutrino as 
summarized in [l3< |. 

Even so, various models have been proposed in which 
singlet "sterile" neutrinos mix in vacuum with active 
neutrinos (i/g, I'p,, Vt)- Such models provide both warm 
and cold dark matter candidates. By virtue of the mix- 
ing with active neutrino species, the sterile neutrinos are 
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FIG. 1: Evolution of various densities as labeled as a function 
of time for a cosmology in which the dark matter begins to 
decay at a redshift = 0.3 (fd « 8 Gyr) with a decay lifetime 
of r = 20 Gyr. The dashed line shows the evolution of rela- 
tivistic particles Pr = p-y + pi- The dot-dashed line shows the 
evolution of total matter pM = Pdm + Pb + ph- The dot-dot- 
dashed line shows the combined evolution of radiation and 
Pbv- The solid line shows the evolution of the total mass- 
energy density . The central dotted line is to aid the eye to 
identify and also the flattening of the total density during the 
decay epoch. This flattening leads to cosmic acceleration and 
apparent dark energy. The present age for this model {t = 11 
Gyr) occurs at (p/pc(now) = 1.0. 



not truly "sterile" and, as a result, can decay. In most 
of these models, however, the sterile neutrinos are pro- 
duced in the very early universe through active neutrino 
scattering-induced de-coherence and have a relatively low 
abundance. As pointed out in [Tsj . however, this produc- 
tion process could be augmented by medium enhance- 
ment stemming from a significant lepton number. Here 
we speculate that a similar medium effect at late times 
might also induce a late time enhancement of the decay 
rate. 

There are several other ways by which such a heavy 
neutrino might be delayed from decaying until the 
present epoch. The possibility of delayed decay by a 
cascade of intermediate decays prior to the final bulk- 
viscosity generating decay was explored in [Tsj . However, 
this possibility was found to be incapable of accounting 
for all of the cosmic acceleration. Here we propose that 
the decay of heavy neutrinos at late times requires one of 
two other possibilities: 1) A late low-temperature cosmic 
phase transition whereby a new ground state causes the 
previously stable dark matter to become stable; or 2) a 
time varying effective mass for either the decaying par- 
ticle or its decay products whereby a new ground state 
appears due to a level crossing at a late epoch. 

The first possibility was considered in Ref . [l^ . A suit- 
able late decaying heavy neutrino could be obtained if the 
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decay is caused by some horizontal interaction (e.g. as 
in the Majoron [23| or familion models). The decay 
within the familion model can be written; 



(24) 



where / is a massless Nambu-Goldstone boson associ- 
ated with a spontaneously broken "family symmetry." 
The lifetime for the heavy neutrino decay then becomes: 



(10^ yr)[100 eV/TO^)3[F/(10^ GeV)] 



(25) 



where F is the scale of the spontaneous symmetry break- 
ing and nil, is the neutrino mass. An appropriately cho- 
sen symmetry breaking scale is required [15j to induce 
the decay at the desired epoch. 

Regarding the second possibility, a number of papers 
have been written ^27] which consider dark matter neu- 
trinos with a time varying mass of (the so-called MaVaN 
model). Those models were proposed ^2T] as a means 
to account for dark energy from self interaction among 
dark matter neutrinos. In our context we would require 
the self interaction of the neutrino to produce a time- 
dependent heavy neutrino mass such that the lifetime for 
decay of an initially unstable long-lived neutrino becomes 
significantly shorter at late times. 

Another possibility might be a more generic long-lived 
dark-matter particle ip whose rest mass increases with 
time 28] . This could be achieved (for example in scalar- 
tensor theories of gravity [1^) by having the rest mass 
derive from the expectation value of a scalar field (p. If 
the potential for (p depends upon the number density of 
ip particles then the mass of the particles could increase 
naturally with the cosmic expansion. This could lead to 
a late-time instability to decay. 

In another proposal [s^l it has been pointed out that 
in string effective theories there can exist a dilaton scalar 
field which couples to gravity matter and radiation. In 
general, different particle masses will have different dila- 
ton couplings. This can lead to dark matter with vari- 
able mass |30|. Moreover, the dilaton couples to radia- 
tion in the form of a variable gauge coupling. This could 
also lead to a variable decay rate of a long-lived dark- 
matter particle such that an initially neglegible decay 
rate quickly accelerates at some late epoch. 

Another possible candidate may be from supersym- 
metric dark matter. It is now popular to presume that 
the initially produced dark matter relic must be a super- 
WIMP in order to produce the correct relic density 0. 
This superWIMP then decays to a lighter stable dark- 
matter particle. One interpretation of a candidate for 
late decaying dark matter here, then might be a decaying 
superWIMP with time-dependent couplings and hence a 
variable lifetime. Alternatively, the light supersymmet- 
ric particle, might also be unstable with a variable decay 
lifetime. It has been proposed [3l| for example that there 
are discrete gauge symmetries (e.g. Zio) which naturally 
protect heavy X gauge particles from decaying into or- 
dinary light particles so that the X particles become a 



candidate for long-lived dark matter. The lifetime of the 
X, however is strongly dependent upon the ratio of the 
cutoff scale (M, « 10^^ GeV) to the mass of the X. 
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1 



10^^ Gyr 



(26) 



From this it is apparent that even a small variation in 
either Mx or could lead to a drastic speed up in the 
decay lifetime. 



IV. SUPERNOVA DISTANCE-REDSHIFT 
ANALYSIS 

Having defined the cosmology of interest we now ex- 
amine the magnitude-redshift relation for type la super- 
novae (SNIa). The apparent brightness of the type la 
supernova standard candle with redshift is given [331 by 
a simple relation for a fiat A = cosmology. The lumi- 
nosity distance in the present cosmological model can be 
written. 



Ho L^o 

+nj,M{.z') + n^,{z') + ^^{z') 



n-Xz') + nx{z') 

1-1/2 



(27) 



where is the current Hubble parameter. The fi^ are 
the energy densities normalized by the present critical 
density, i.e. ^i{z) = 87rGpi{z)/3HQ. ilh is the closure 
contribution from the decaying heavy cold dark-matter 
particles. Their decay is taken here to produce light neu- 
trinos or other noninteracting relativistic particles 
while Q-y is the contribution from normal relativistic mat- 
ter. Note that f2hj and l^i each have a nontrivial 
redshift dependence due to particle decays, while flj{z) 
varies as (1 -I- z)~* and stable dark matter and baryons 
fl-DM{z) + rib(z) obey the usual (1 -I- z)~^ scaling. 

Figures O and [3] compare various cosmological mod- 
els with some of the combined data from the High- 
Z Supernova Search Team and the Supernova Cosmol- 
ogy Project [H, [3^. These figures shows the evolution 
of the relative distance modulus as given in the usual 
way A(m - M) = 5log[DL/ DL{flk = 1)], where the K- 
corrected magnitudes m = M + 5 log Dl + 25 are plotted 



relative to to a fiducial ilk 
mology, for which 



l/{aoHo)- — 1 open cos- 



1) 



cjl + z) 
2H„ 



z + l 



(^ + 1) 



(28) 



As a measure of the quality of the different models. 
Table 1 summarizes the relevant model parameters and 
reduced goodness of fit. The quantity in column 
4 is defined as the present sum of nonrelativistic matter. 



M 



nh{z = 0) + nuuiz = 0) + f7b(z = 



(29) 




r (Gyr) 


Zd 


C 






xi 


Age (Gyr) 


> 20 


0.30±0.04 


1. 


0.50 ±0.03 


0. 


1.32 


11.0 ± 0.1 


(ACDM) 




1. 


0.31 


0.69 


1.14 


13.0 


(CDM) 




1. 


0. 


1. 


3.23 


9.2 


20 


(no delay) 


1. 


0.16 


0. 


3.93 


8.7 


20 


0.3 


1.6 


0.30 


0. 


1.13 


12.0 


5 


0.26 


4.0 


0.27 


0. 


1.13 


11.9 



TABLE L Parameter sets for various fits to the SNIa 
luminosity-redshift relation for Hq = 71 km s~^ Mpc~^ and 
fib = 0.044. In the decaying (finite r) models no stable dark 
matter was assumed (i.e. SIdm = 0). 



FIG. 2: Fit to the SNIa magnitude-redshift relation for an un- 
enhanced bulk viscosity from late decaying dark matter (solid 
line) compared with the observations [3^ (points). The fit 
corresponds to Zd = 0.3 {td ~ 8 Gyr), and r > 20 Gyr with 
an age of 11 Gyr. Data and lines are plotted relative to a 
fiducial Q.k = 1 open cosmology. For comparison, the dashed 
line is for a flat Q,m = ^ matter-dominated cosmology. The 
dot-dashed line shows a flat standard ACDM cosmology with 
JIa = 0.7 and JIm = 0.3. The Dash-double-dot line shows a 
model as in [l^ with undelayed particle decay. 



Zj = 0.3 ; T = 20 


C= 1.6 




0.0 




0.7 


■ no delay 




Zj = 0.24 ; T = 5 


C = 4.0 




FIG. 3: Same as Figure [J] but for additional parameter C — 
1+aTe due to possible higher order effects in the bulk viscosity 
(. The solid line corresponds to Zd ~ 0.3 {td ~ 9 Gyr), and 
C = 1.6 with r = 20 Gyr and an age of 12 Gyr. The dash- 
dash-dotted line corresponds to Zd = 0.26 (td ~ 9 Gyr), and 
C = 4.0 with r = 5 Gyr and an age of 12 Gyr. 



The solid line on Figure [2] shows a fit to the SNIa data 
for a model with a fixed bulk viscosity coefficient as given 
in Eq. [211 The solid line and dash-dash-dotted lines on 
Figure [3] show fits to the data when an enhanced C is 
allowed due to possible higher order terms in the trans- 
port equation. The dotted line on both figures shows the 
fiducial Ofc = 1 open cosmology. For comparison, the 



dashed line on both figures is for a fiat f^M = 1 matter- 
dominated cosmology, while the dot-dashed line shows 
a flat standard A-|-cold dark matter (ACDM) cosmology 
with = 0.7 and f^M = 0.3. 

The fit with the bare bulk viscosity coefficient (Eq. YH^ 
on Figure m corresponds to Zd = 0.3 {td ~ 8 Gyr), and 
T > 20 Gyr with an age of 11 Gyr. The Dash-double- 
dot line shows a model as in with undelayed particle 
decay. As discussed in a model without delayed 

decay actually produces a worse fit than even a matter- 
dominated model. The reason is simply that the forma- 
tion of a radiation-dominated universe by particle decay 
causes the energy density to diminish even faster with 
scale factor than a matter dominated cosmology, even 
with the help of the bulk-viscosity term. On the other 
hand, by delaying the onset of particle decay and the as- 
sociated bulk viscosity until late, it is possible to produce 
some of the features of a ACDM model without invoking 
a A. 

As a practical matter, it turns out that with the un- 
enhanced bulk viscosity in Eq. the best fit requires 
the largest bulk viscosity. From Eq. ([H]), however, it 
is clear that large ^ occurs for large Tg in the limit of 
T >> (3i/)~^. On the other hand, an infinite lifetime 
should imply no viscosity as there is no decay. We resolve 
this dilemma by imposing an arbitrary cutoff in the decay 
lifetime of a couple of Hubble times (~ 20 Gyr). There is 
not much change in the goodness of fit for larger times. 
The decay lifetime then is not a parameter, it is simply 
large (> 20 Gyr). The only parameter in this model is 
therefore the time (or redshift) at which the decay begins. 

On the other hand, when the enhancement factor C 
is introduced as shown of Figure [31 equivalent fits can 
be obtained with almost any value for r. So again, r is 
not a parameter, or more precisely, it has a degeneracy 
with C, such that any decrease in r can be offset by an 
increase in C. Moreover, because the total bulk viscosity 
can be increased, the goodness of the fit is substantially 
improved, even after allowing for the introduction of an 
additional degree of freedom. 

One interesting feature of the bulk- viscosity models ap- 
parent in both Figures [H and [31 is that the magnitude- 
redshift relation decreases more rapidly (brighter appar- 
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ent magnitude) than in a standard ACDM cosmology. 
This can be traced to the higher matter content in the 
past which leads to a more rapid deceleration during the 
matter-dominated epoch. In principle, this feature could 
allow one to distinguish between the two cosmologies as 
more SNIa data is accumulated at highr redshift. For 
now, the data slightly favor the bulk-viscosity models. 

Of particular interest regarding Figures[21[31 and Table 
1 is the fact that a A = model can produce a reduced 
which is as good as the best fit standard ACDM cosmol- 
ogy. This is based, however, upon a parameterization 
of the effects from higher-order terms in the transport 
equation, and needs to be verified. Even so, it is at least 
possible that this late-decaying model for bulk viscosity 
realizes the present cosmic acceleration without the need 
for a cosmological constant, and hence, is a viable alter- 
native to the ACDM model based upon fits to the SNIa 
data. 

The cosmic age and the growth of large-scale structure 
place additional constraints on a late-decaying cosmology 
[l5|. However, for the optimum fits derived here the on- 
set of the decay occurs close enough to the time when a 
normal ACDM cosmology begins to become A dominated 
that there is not much difference in the implied cosmic 
age. The ages deduced here are lower ( 11-12 Gyi) than 
the age (13.8!^g'2 Gyr) deduced from an analysis [3^| of 
WMAP data based upon a ACDM model. However, a 
direct comparison with the WMAP age is not meaning- 
ful unless the CMB analysis is redone using the present 
cosmological model. This we plan to do in a future work. 
Nevertheless, the present model is consistent with the in- 
dependent constraints from the age from the oldest stars 
(13±2 [11]) and globular clusters (11±3 (Hi). Moreover, 
even if the age is low, this problem could be alleviated in 
open cosmological models with bulk viscosity [1] . 

The formation of large scale structure in the present 
model, however, might be a more serious problem. Even 
though the expansion rate is not much different than in 
the best fit ACDM model, in this model the dark matter 
content is higher in the past. This may lead to excess 
of large scale structure at early times. Also, the current 
matter content in some of the fits, flM = 0.27 — 0.50, 
is high compared to the value derived from the WMAP 
[S^ l analysis VIm = 0.26to o3- Both of these issues could 
be resolved, however, by considering an open cosmology. 
Clearly, this is something which needs to be investigated 
and in a future work we will examine these issues. We do 
note, however, another positive feature of cosmic struc- 
ture in these models [ISJ. The formation of large scale 
structure in a cosmology with decaying dark matter can 
lead to a flattening of the dark matter density profiles 



consistent with observations [3' 



V. CONCLUSION 

We have considered models in which the present cos- 
mic acceleration derives from the temporary insertion of 
dissipative mass energy due to the bulk viscosity created 
by the recent decay of a cold dark-matter particle into 
light (undetectable) relativistic species. As illustrative 
examples we have considered initially matter dominated 
flat (A = 0) cosmologies plus late-time particle decay. 
We find that models with bulk viscosity from late-time 
dark-matter decay are consistent with observations of the 
supernova magnitude-redshift relation, and ages from the 
oldest stars and globular clusters. We argue that it will 
likely satisfy other constraints as well. 

Moreover, there is a difference in the SNIa magnitude- 
redshift relation for this cosmology compared to the stan- 
dard ACDM model. This is because the deceleration is 
faster at high redshift due to a higher matter content 
during the matter-dominated epoch. Thus, as more data 
are accumulated at the highest redshifts, it may be possi- 
ble to distinguish between this cosmology and a standard 
ACDM model. For now, however, there is sufficient suc- 
cess in the present model to motivate further work. In a 
subsequent paper we plan to consider higher order terms 
in the transport equation in detail as well as the effects 
of such late decays on the observed power spectrum of 
the cosmic microwave background [34| . and the growth 
of large scale structure. 

Ultimately, of course, one must decide whether the 
dilemma of a cosmological constant with all of its difficul- 
ties is less palatable than the dilemma of a bulk viscosity 
produced by the delayed onset of dark-matter decay. For 
now, however, our purpose has simply been to establish 
that such a possibility exists and that it warrants further 
investigation. 
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